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Abstract 

We present a method to construct a symplecticity preserving renor- 
malization group map of a chain of weakly nonlinear symplectic maps 
and obtain a general reduced symplectic map describing its long-time 
behaviour. It is found that the modulational instability in the reduced 
map triggers random wandering of orbits around some homoclinic-like 
manifolds, which is understood as the Bernoulli shifts. 
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1 Introduction 



Hamiltonian systems with more than one degree of freedom can exhibit 
trajectories with complex behavior. Discrete symplectic maps are of inter- 
est since the Poincare surface of section leads naturally to a mapping of the 
dynamical trajectory onto a subspace of the phase space of the Hamilto- 
nian flow. Symplectic maps allow much easier numerical calculations of the 
motion than Hamiltonian flows, which are particularly convenient for study- 
ing systems of many degrees of freedom. In a two-dimensional symplectic 
mapping, which corresponds to Poincare mapping of continuous-time Hamil- 
tonian systems with two degrees of freedom, the main ingredients of this 
chaotic motion are hyperbolic type orbits and their invariant manifolds. The 
existence of transversal intersection between stable and unstable manifolds 
leads to non integrability of Hamiltonian systems. The area of the lobe en- 
closed by stable and unstable manifolds represents the flux from inside (resp. 
outside) the separatrix to outside (resp. inside). The manifolds which have 
codimension one separate the space into disjoint regions. 

For symplectic mappings of four or more dimensions, there are some stud- 
ies about the homoclinic bifurcation near fully hyperbolic fixed points |3J 
Although the dimension of unstable and stable manifolds of a fully 
hyperbolic fixed point takes a maximum possible value, the transversal in- 
tersection occurs on a nearly one-dimensional homoclinic submanifold and 
exponentially small splitting of the homoclinic manifold leads to very weak 
chaos, which is almost invisible through naive numerical experiments. 

A chain of weakly coupled nonlinear oscillators may present a model for 
dynamics near a fully elliptic fixed point. The equilibrium state of almost 
independent oscillators corresponds to an fully elliptic fixed point in the phase 
space and periodic orbits or tori encircle the elliptic fixed point. Resonant 
tori may be disrupted by the well-known resonance mechanism. If a torus is 
hyperbolic, unstable and stable manifolds attach to the torus as a "whisker" 
H. Some dynamics of such whiskers attached to tori is studied near the 
resonance junction by means of the naive averaging method |J . 

In the case of a chain of weakly coupled nonlinear oscillators, where both 
coupling and nonlinearity are weak, tori may be destabilized by another 
mechanism called the modulational instability. Here, we study about such 
modulationally unstable tori and their unstable manifolds. We start with 
a chain of symplectic maps, each of which represents a weakly nonlinear 
oscillator. In order to avoid special dependence of results on the specific 
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form of nonlinearity, our investigation is focused on long-time asymptotic 
behavior of dynamics of the system. However, the conventional asymptotic 
methods such as the averaging method or the method of the multiple-time 
scales may not be immediately applicable to discrete systems. Instead, we 
employ the perturbative renormalization group (RG) method |§0, which 
can apply to a discrete system and leads to a reduced map as the RG map 
||. However, a naive reduced map does not preserve symplectic symme- 
try and fails to describe long time behavior of the original symplectic map. 
For two-dimensional maps or special four- dimensional maps, the symplectic- 
ity is recovered through a simple "regularization" procedure ||. For higher 
dimensional symplectic maps, we do not have a general procedure of "regu- 
larization" . The first purpose of this paper is to develop a general procedure 
of "regularization" of the RG map for a higher dimensional symplectic map. 
Our regularization procedure is based on the symplectic integration method 
for a continuous Hamiltonian flow and includes the previous regularization 
procedure as a special case. By means of this regularization procedure, we 
can derive a reduced symplectic map as the RG map, which asymptotically 
approximates a chain of weakly nonlinear symplectic maps. 

The reduced map obtained here is general in the similar sense that the 
nonlinear Schrodinger equation asymptotically describes a general weakly 
nonlinear dispersive wave. In fact, the reduced map is found to be a space- 
time discrete version of nonlinear Schrodinger equation. Since the reduced 
map is not integrable, while its continuous limit is integrable, it is worth 
studying chaotic solutions of the reduced map, which is the second purpose 
of this paper. The reduced map has a periodic solution, which is destabi- 
lized by modulational perturbation under some conditions (the modulational 
instability). Then, the periodic solution becomes hyperbolic and is associ- 
ated with unstable and stable manifolds, whose transversal intersection is the 
origin of chaos (the homoclinic chaos). It should be mentioned that various 
discrete versions of the nonlinear Schrodinger equation are studied in connec- 
tion with numerically induced chaos of the nonlinear Schrodinger equation 
pifl[]n] • There, the systems of a large degrees of freedom are investigated 
in connection with numerical integration of the nonlinear Schrodinger equa- 
tion and only few attentions are directed to behavior of chaotic solutions 
themselves. 

Here, we are mainly interesting in chaotic solutions near the onset of 
symplectic chaos, where the phase space of the map has a dimension higher 
than two. Near the onset of the homoclinic chaos, an unstable manifold is 
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found to be very close to a stable manifold and a homoclinic-like structure 
is approximately preserved. In a two-dimensional symplectic map, the ho- 
moclinic chaos close to the onset is hardly seen since a splitting distance 
between unstable and stable manifolds is small beyond all order of the per- 
turbation parameter. However, in a higher dimensional map considered in 
this paper, the weak chaos is shown to be visualized as a random rotation of 
some homoclinic-like orbits, which is characterized by the Bernoulli shift. 

In the section 2, some reduced maps are derived from a chain of symplec- 
tic maps as the regularized RG map. The conditions of the modullational 
instability in the reduced maps are presented and unstable manifolds of un- 
stable periodic orbits are depicted in the cases of two, three and four sites in 
the section 3. In the section 4, a random rotation of homoclinic-like orbits 
is observed for three- and four- site maps, which is interpreted as the the 
Bernoulli shift of three or four numbers. 

2 Symplectic Map Chain and Reduced Map 

In this section, we consider a symplectic map chain and present a method to 
construct a symplecticity preserving renormalization group map ( a reduced 
symplectic map) from the symplectic map chain. 

2.1 Symplecticity Preserving Renormalization Group 
method 

We give a method to derive a symplecticity preserving reduced map from a 
nonlinear symplectic map of three or more sites, which is non-integrable even 
in the time-continuous limit. 

4 +1 = X]+T P ] +1 , (1) 

PT 1 = p] + r(-n 2 x] + e{v(x] +1 -2x] + x]_ 1 )-a(x]) 3 }), (2) 

with N(> 3)-periodic boundary condition : (x™ +A r,p™ +A r) = Here 
x™ and p™ are real canonical variables on site j at time n, while Q, a, u, and 
t are real parameters and e is a small perturbation parameter, < £ < 1. 
Eliminating p -variables from Eqs. (|1|) and (fj), we obtain a map for x- 
variables: 

x] +1 - 2 cos(6)x] - x^ 1 = ST 2 {v(x n j+l - 2x] + x]^) - a(x™) 3 }, (3) 
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where cos(#) = (1 — Q 2 r 2 /2). Let us construct a reduced map for < e <C 1 
by means of the RG method. Expanding a solution of Eq. (0) as x™ = 

xf 0) + exf 1] + 0(e 2 ), j = 1, 2, .., we have 



x 



n(0) 
3 



Ajexp(—i8n) + c.c, 



Cxf 1] = t 2 \v A 2 Aj - 3a 1^1%- 1 exp(-i6n) + c.c, 



where 



Cxf l) = X ; +1{1) -2cos6-xf 1) +x;- 1{1 \ (4) 
A 2 A] = A] +1 -2A] + A]_ V (5) 

and c.c. denotes complex conjugate to the preceding terms; Aj is an integra- 
tion constant associated with a site j. The first order solution has a secular 
term with respect to n. 



n(i) ir 2 

rr> ■ r iy-\ 

i ~ 2sin0 



jz/ A 2 Aj — 3a\Aj\ 2 Aj^ exp(—i9n) + c.c. + n.r., 



and n.r. denotes non-resonant terms. In order to remove this secular term, 
we define a renormalization transformation Aj i— > A" by 

^ = A i + ^{n^ (v A 2 A, - 3a|4-|%) } + 0(e 2 ). (6) 
A discrete version of the RG equation is constructed by taking difference of 

An 

Substituting the expression for Aj in terms of AJ defined by Eq. (|6|) into 
Eq. (0), we can eliminate the secular term up to 0(e) and obtain a naive 
RG map 

A T = A l + £ {^§( v A ? A ^ ~ ^ A ^ A i) }> ( g ) 

or 

/in+l An ■ 

This naive RG map does not preserve symplectic symmetry. To recover 
symplecticity of Eq. (M), we apply the symplectic integration method to the 
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continuous-time limit (r — > 0) of the naive RG map (Q). The continuous- 
time limit takes the following symplectic form, which is also derived from the 
continuous-time limit of Eqs. ([[]) and (Q) in the Appendix A, 

dA, If a2 tn ,„„ l2 „„\ dH RG dA* dH RG . , 

sin(#) Rr 
Here, 12 = as r — > and if is given by 

T 



H RG = + 



•J I ' 



2ft ^ 1 3+1 

3 

An Jl 

i 

The Hamiltonian flow defind by if a can be solved as 



Aj(t + r') = Aj(t) exp ( - -i^|^-(i)|V 



where r' is an arbitrary constant and we have 



AfW^ xp (-^|^|V). (ll) 



The flow defined by H v can be discretized by the symplectic implicit midpoint 
rule [0: a general Hamiltonian flow 



is discretized as 
In our case, the scheme gives 

^r 1 = ^ n + 4^ ( a ? ^ n+1 + A ?^ n ) > (12) 

or a temporally explicit form 



Combining two symplectic transformations (JTTJ) and ([T3|) , we get a sym- 
plectic scheme for H = H v + H a as 



For brevity, we rewrite the above symplectic RG map as 

(l - iT A 2 AA] +1 =(l + iT A 2 ) exp (iQ\A]\ 2> \A^ (14) 

where 

z/t 2 -3«r 2 „ 

4 sin & 2 sin d 



Which has a conserved quantity ^ ■ |A™| 2 (see Ref. ||13|| ). Both temporally 
and spatially explicit form of Eq. flTl]) is possible. For instance, the three-cite 
map is written as 

_ 1 - 3iT 
~ 1 + 9T 2 









)l 









where 



^exp^QI^I 2 )^ 



An approximate solution of the original map ([!]) and @ is expressed in terms 
of the renormalized amplitude as x™ ~ A"exp(— iQn) + c.c. If the coupling 
parameter v vanishes, the present regularization process is nothing but a 
simple exponentiation procedure introduced in ||. Note that the Eq.(|T4]) 
with N 3> 1 is identical to a difference scheme for the nonlinear Schrodinger 
equation [13|]. In fact, the time-continuous limit of the RG map ( |T4T) is just a 
spatially discretized nonlinear Schrodinger equation (see Appendix A). The 
map QHD also admits some oscillating solutions including a spatially uniform 
one. 



2.2 Nonlinear symplectic map with two sites 

Let us consider a nonlinear symplectic map with two sites. This system is 
separately discussed since the continuous-time limit (r — > 0) of its RG map 
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is integrable. 

x^ 1 = x" + rp™ +1 , =x" + rp™ +1 , (15) 

p y+i = p? + r{-fi 2 a:? + e(i/(x2-xi)-az?)}, (16) 

p«+i = p^ + r{-fl^ + e(^i-x 2 )-a^)}. (17) 

We rewrite Eqs. (|T5|) - (|17|) as follows, 

-2cos^-^ + x^ 1 = er 2 |z/(x™ - x") -q^) 3 }, 

-2cos^-^ + x™- 1 = er 2 {z/(£™ -x") -«04) 3 }, 

where 6 1 is defined in Eq. ([3]). Setting x 1 - ~ A™ exp(— iOn) +c.c. and following 
the same procedure as in the case of the three- or more- site model, we reach 
a regularized RG map up to 0(e) 

(1 - iTLj)AJ +1 = (1 + iTLj) exp(tQ\A]\)A], (j = 1, 2) (18) 

where an operator Lj is defined by 

L\A\ = -A\ + A$, L 2 A% = A™ — A". (19) 

The detail derivation of this regularized RG map is shown in Appendix B. 
It is worth noting that the map ([Tj^) has not only symplectic symmetry but 
also a conserved quantity 

£ l^l 2 = £ KT- ( 2 °) 

j=l,2 j=l,2 

The map (|18D has oscillating solutions 

Ai = A 2 = A°exp(iQ\A°\ 2 nj, (21) 
or 

A" = -A 2 = A°exp (iQ|A°| 2 ™), (22) 

where A is a constant. The family of these solutions lies on a one-dimensional 
torus in the phase space that spanned by A™ and A" *. The (modulational) 
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stability of solutions (|2"I|) and (^) is of interest because the stability de- 
termines hyperbolicity or ellipticity of the torus. If the torus is hyperbolic, 
unstable and stable manifolds attach to the torus like "whiskers" 0] and we 
call such unstable and stable manifolds whiskers in this paper. 

Behavior of such whiskers is the main object to investigate in the following 
sections 3 and 4. Let us show effectiveness of the regularized RG map by 
numerical calculations. In Fig. [3] , trajectories obtained from the naive RG 
map (H) and the regularized RG map (|14D are depicted to be compared to 
an exact trajectory of the original map (Eqs. (HD-(||])). All initial points are 
set very close to a hyperbolic torus. The trajectory given by the regularized 
RG map is in good agreement with the exact solution. 

Owing to the conserved quantity Q2l]|), the time-continuous limit (r — > 0) 
of the RG map (|i~8| ) becomes integrable and analytical expression of whiskers, 
which is exactly homoclinic in this case, is given in the Appendix C. 



3 Whiskered Tori 

In this section, we analyze the modulational instability of spatially uniform 
oscillating solutions (tori) of the reduced maps derived in the previous section 
and introduce whiskered tori, that is, the hyperbolic structure near periodic 
manifolds. It is found that whiskers retain the homoclinic structure approx- 
imately near the onset of the modulational instability, where the homoclinic 
chaos is weak. 

3.1 Two-site case 

Let us analyze the reduced symplectic map with two sites, whose continuous- 
time limit is integrable. 

One of tori of Eq. ( |T8| ) is just a spatially uniform solution 

A] = A°exp(tQ\A°\ 2 n), (j = l,2). 

To study the phase space structure around the solution, we consider a per- 
turbed solution of the form 

A] = A°exp(zQ|A°| 2 n)(l+^ n ), (23) 
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where |/i"| 2 <C 1. Substitution of this expression into and retaining first 
order terms in /x™ , yield 

1 + iT —iT \ ( fi n+1 \ 2 ( 1 — iT iT \ ( u n 

-it i + iT = ( 1 + i Q\ A °\ 2 )y lT 1 _ iT jy^n 

1 2 



and complex conjugate to the above expression. The eigenvalues of Eqs.( 
are 

1, 1, \±=(3±y/P*-l, (25) 

where (3 = (1 - 4T 2 + 4Q\A°\ 2 T) / (1 + AT 2 ). We obtain conditions for the 
instability 

\A°\ 2 > ^, (26) 

or 

I-4T < (27) 

which comes from (3 2 > 1. Hereafter, we concentrate on the condition (PB|), 
under which the torus becomes hyperbolic for the amplitude larger 
than a critical value (y/2T/Q). From the spectrum of eigenvalues (p5|) , we 
find that whiskers consists of one-dimensional unstable manifold (departing 
whisker) and one-dimensional stable manifold (arriving whisker). To con- 
struct whiskers numerically, we transform A n - to a™(j = 1,2) where 

A] = a]exp{iQ\A°\ 2 n). 

The torus with departing whiskers are depicted in Figs 0. Even if \A°\ is 
considerably larger than the critical value, whiskers are found to be almost 
homoclinic and their behavior appears to be regular. This apparent regularity 
reflects the integrability of continuous-time limit of the map with two sites 
and is sharp contrast to the case of three or more sites as studied in the 
subsequent sections. Since difference between the map (|T8| ) and its integrable 
continuous-time limit is very small (beyond all order), the weak homoclinic 
chaos may be almost invisible in the present range of numerical calculations. 
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3.2 Three-site case 



Let us consider the regularized RG map with three sites When the 

site number is odd, general stability analysis is possible and is provided 
in Appendix D. Similar linearized analysis around a uniformly oscillation 
solution as in the previous subsection leads to instability conditions: 

|AT > ^, (28) 

or 

i- 40 ! 2 < m (29) 

As well as the two-site model, we consider the condition (|28|) only The 
spectrum of the linearized map is given as 



A a = (3 a ± ^(31 - 1, (a = 0,1) (30) 

Po = 1, (singlet) 

1 + 6Q|A°| 2 T-9T 2 
ft = 1 + QT 2 , (doublet), 

The spectrum ( j30|) shows that the dimension of both departing and arriv- 
ing whiskers are two for three-site case. The tori with its unstable manifolds 
are depicted in Figs. These figures show that the whiskers' behavior is 
different from that of two-site case. The difference may be attributed to one 
of integrability property of their continuous-time limit. The greater is the 
amplitude of tori, the more complicate is the whisker's behavior. When the 
amplitude is just above a critical value (y/ST/Q) for the instability, whiskers 
are found to be still homoclinic-like as shown in Fig [5] (a) . 



3.3 Four-site case 

We analyse the modulational instability of a uniformly oscillating solution of 
the regularized RG map with four sites (fTJ) separately since general analysis 
is possible only for the odd-site case (Appendix D). 

The linearized map around a uniformly oscillating solution (a torus) has 
the following eight eigenvalues 

A a± = /3 a ± V/^ 2 - 1, (a = 0,1, 2) 
Po = 1, (singlet) 
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1 + 4T2 , (doublet) 

l + 8g|A°| 2 T-16T 2 
& = 1 + 16T2 , (singlet). 



The eigenvalues give the following instability conditions of the torus. 

2T 47 
(doublet), |A°| 2 >^ 



IT AT 
> (doublet), \A Q \ 2 > (singlet), (31) 



or 

\ A "\ 2 <-^f (doublet), |A°| 2 <-^ (singlet), (32) 

where we consider the conditions fl3~ID only as before. The spectrum of eigen- 
values gives the dimension of whiskers associated with the torus: 

2T 

• |A 1 2 < — • • • no whiskers (elliptic), 

2T AT 

< | A 1 2 < — — ■ • ■ 2 dimensional departing and 2 dim arriving 



Q Q 

AT 

~Q 



whiskers and ellipticity is equal to 2 dimensions, 

AT 

\A°\ 2 > — • • • 3 dimensional departing and 3 dimensional 

arriving whiskers. 



The torus with its unstable manifold is depicted in Fig. [|. As similar to 
the three-site case, the greater is the amplitude of torus, the more complicate 
is whiskers' behavior. 



4 Random Rotation of Homoclinic-like Or- 
bits 

For slightly above the critical amplitude, we study chaotic orbit's behavior 
in the three- and four-site cases. In such cases, an orbit wanders randomly 
around some homoclinic-like manifolds (whiskers) and it looks like a random 
sequence of homoclinic-like orbits as shown in Fig. [5]. Each homoclinic- 
like orbit only differs in which site takes the largest amplitude in the orbit. 
In other word, a random rotation of site-numbers occurs in a sequence of 
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the same homoclinic-like structure. This random rotation of site-numbers 
is understood by means of the Bernoulli shift. Let us characterize each 
homoclinic-like orbit by the site-number of the largest amplitude so that an 
orbit is represented by a sequence of site-numbers 

10210120- ■■ (resp. 112301 ■••), etc. 

corresponding to the sequence in Figs. |6| and [7], where indicates the site 
number 3 or 4. This random sequence of numbers is conceived by the 
Bernoulli sequence, which is generated by a simple map 

w^ 1 = 3w n (resp. w n+1 = Aw 11 ), mod 1. 

A periodic sequence of homoclinic-like orbits corresponds to a periodic so- 
lution of the Bernoulli map for a rational initial value on the interval [0, 1), 
while a random sequence is generated for an irrational initial value. 

Now, we present some number sequences of the homoclinic-like whiskers 
obtained by numerical iterations of the three-site case and the four-site case 
respectively in Fig. [8[ These figures suggest that the sequences of homoclinic- 
like orbits is not periodic but random. It should be noted that such a random 
rotation of homoclinic-like orbits does not appear in the two-site case. 

5 Conclusion 

We present a regularization procedure to preserve the symplectic structure 
of the RG map near a fully elliptic fixed point of a chain of weakly nonlinear 
symplectic maps. The regularization is accomplished by extended exponen- 
tiation of the naive RG map and we derive a general reduced symplectic map 
as an asymptotic RG map. 

Analyzing the modulational instability of a uniformly oscillation solution 
of the reduced map, we find a hyperbolic torus with whiskers. We observe 
that whiskers retain a homoclinic-like structure in the case the amplitude of 
the uniformly oscillating solution exceeds slightly the critical value. For three 
or more-site cases, we find a random sequence of the homoclinic-like orbits. 
The orbit is symbolized by a random rotation of site-numbers generated by 
the Bernoulli shifts. Although this chaos may be produced by the homo- 
clinic mechanism and is very weak, it is easily visible as a random sequence 
of homoclinic-like orbits. Such easily visible irregularity of homoclinic-like 
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orbits may be general in high- dimensional chaos and give an important vis- 
ible tool to understand the onset of chaotic dynamics in high-dimensional 
spaces. 



Acknowledgement 

The present work is, in part, supported by the Japan Society for Promotion 
of Science, Grand-in-Aid for Scientific Research (C) 13640402. 



Appendix 

Appendix A: Continuous- time RG equation with three 
or more sites 

Here, we give an RG equation of continuous-time limit (r — > 0) of the 
weakly nonlinear chain (|l|) and (|2|), which is a spatially discretized nonlinear 
Schrodinger equation. Continuous-time limit of Eqs. ([!]) and (Q) reads 

dxj dH 

dpj r^i { ( \ o "| dH 

— = -fi x j + E^v[x j+1 -2x j + x j - 1 ) -a Xj j = - — , 



i i 

Expanding Xj as Xj = xf^ + sXj + 0(e 2 ), a naive perturbed solution is 
obtained as 



x 



= Aj exp(-iQt) + c.c. 



xy = —[uA^Aj- 3a\Aj\ 2 Aj) exp(—i£H) + c.c. + n.r 



where Aj is a complex valued integration constant for j-th site. To remove 
secular terms proportional to t, we define the renormalization transformation 
Aj i— > Aj (t) by collecting all terms proportional to the fundamental harmonic 
exp(— iflt) 

Aj(t) = Aj + (y A 2 A, - 3a|^|%) + 0(e 2 ), 
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From this transformation, the prescription of the RG method gives an RG 
equation 



dAj _ . 1 



(vA'Aj-dalAjl'Ajy (33) 



This RG equation automatically retains the symplectic property and is known 
as a spatially discretised nonlinear Schrodinger equation, which is also non- 
integrable except the case of two sites. 

Appendix B: Derivation of a reduced map of two-site 
nonlinear symplectic map 

In this subsection, we derive a regularized RG map from a nonlinear two-cite 
symplectic map chain. 

Expanding and x 2 as x™ = x™^ + £%j + 0(e 2 ), (j = 1, 2), we have 
the following naive perturbative solution of Eqs. (|l^) - (JTj 



x 



"(0) 

j 

n(l) 



Aj exp(— i9n) + c.c, 



x 



x 2 



t 2 ^v(A 2 - Ai) - 3«|Ai| 2 Ai| exp(-i9n) + c.c. + n.r., 
t 2 ^v(Ai - A 2 ) - 3a\A 2 \ 2 A 2 ^ exp(-i0n) + c.c. + n.r. 



The renormalization transformation Aj i— > A" to remove the secular terms 
in the coefficient of exp(— iOn) is given by 

A? = Ai + eu^—{v{A 2 - A t ) - SalA^AA + 0(s 2 ) 
2 sin p I J 

A" = A 2 + en^—lv^ - A 2 ) - 3a|A 2 | 2 A 2 ) + 0{e 2 ), (34) 
2 sin ok ) 

from which we have a naive RG map up to 0{e) 

A] = A 1 +e 1 ^(vL J A 3 -3a\A]\ 2 A]), (j = l,2) (35) 

where, LjAj is defined by Eq. (|i~9|). The naive RG map fl35| ) should be 
regularized by means the symplectic integration method. Then, we get a 
regularized RG map: 

(l - iTL^A n + l = (l + iTLjj exp (iQ\A^A h (36) 
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which has the following matrix form useful to numerical calculations: 

( 1 + iT -iT \ ( A™ +1 \ ( 1 - iT +iT \ f B[ l 
\ -iT 1 + iT J \ A T 2 l+l ) ~ \ +iT 1 - iT ) \ Bl 

By inverting the matrix in Eq. fl3"T|) , we have 



.(37) 



A™ +1 \ 1 - 2iT ( 1 2iT \ ( B" 



A™ +L J l + 4T 2 V 2iT 1 J \ B 2 

Appendix C: Continuous- time RG equation with two 
sites 

In order to picture a homoclinic manifold resulting from the modulational 
instability clearly, let us analyze a continuous-time RG equation with two 
sites, which is integrable. In this case, the continuous-time RG equation has 
the following form (see Eq.([33|)) 

dAi is ( ~ ~ . ~ ,,~ I 



(IA2 is c 



- A 2 ) -3a\A 2 \ 2 A 2 y 



Since the regime where the modulational instability occurs is of interest, we 
suppose v = — \u\ < and a > and introduce new variables: 

/ , ! £ ^ /3a ( ivet ~\ ~ , . „ „. 



Then, we have 

dQ 



d / = -i(Q 2+ m*Qi) 

dQ2 -*(Qi + |Q 2 | 2 Q 2 ) 



dU 


dQ\ 


&H 


dQl' 


dt' 


dQi 


&H 


dQ\ 


on 


dQl 1 


dt' 


dQ 2 



dt' 

n = -*{q\Q2 + QiQl + (IQ1I 4 + |Q 2 | 4 )/2}, 

which has a uniformly oscillating solution 

Q 3 (t) =Q (0)exp{ -z(l+ |Qo(0)| 2 ) t) (j = 1,2). 
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To construct a homoclinic manifold of the unstable periodic solution, we 
transform Qj into 

Q j (t')=a j (t)e^{-i(l + \Qo(0)\ 2 ) *'}. 



and have 



w 

da 2 
~dt 



da 2 . r f 



z[{|a!| 2 - (1 + |Qo| 2 )" 


>ai + a 2 


z[{|a 2 | 2 -(l + |Q H 


>a 2 + ai 



(38) 
(39) 



This system has two conserved quantities, i.e. the Hamiltonian and |cti| 2 + 
| a-2 1 2 , and is integrable. After some manipulations of Eqs.(^) and (|39"D, we 
obtain an explicit expression of homoclinic manifold: 



a x = bVl + Texp{i(Q + A)/2}, 
a 2 = by/ 1 - rexp{i(0 - A)/2}, 



and 

m 

A(0 

e(0 



-V^l sech(-2v / ^ r Tt / ) 



Cos" 



lf l-b 2 T 2 /2 
1 - b 2 /2 B 



b 2 -l y/l^B 2 



Tan 1 



1 - B 2 



-tanh (2V& 2 - It') ^ + const. 



where 5 = ly/W^ljb 2 and 6 = |Q (0)| > 1. 



Appendix D: Modulational instability for the case of 
odd sites 

In this section, we provide instability conditions for a uniformly oscillating 
solution of the reduced map ( |HD with odd sites by means of the Fourier 
analysis |13 

set as 



We substitute Eq. (|23|) into Eq. (0), where perturbation is 



2tt 



(40) 
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r N-l N -1 JV-1 N - li 

aG |-^'-^ + 1 '---'°'---'^ }' 

and N is an odd number (N > 3). Assuming the periodic boundary condition 
/•fj+jv = A*", we have a linearized equation for /x™. 

Af 1 \, _ ( d a (i+tQ\Ar) d al Q\Ar \f pz \ uu 
d-\-iQ\Ar) c(i-^T)JUr ( 1 

where a ^ and 

_ 1 -4iTsin 2 (7ra/A^) 
a = 1 + UT sin 2 (na/N)' 

The eigenvalues of Eq. (f£T|) are 

A a± = ft ± V^ 3 !, (42) 

where 

1 + 8QT|A o | 2 sin 2 (7ra/A0 - Ut sin 2 (na/N) 

(3 a = — G R. 

1 + (4Tsin 2 (7ra/iV)J 

Note that A a+ A a _ = 1. The uniformly oscillating solution is destabilized 
when 

Pa > 1 or /3 a < 1, 

that is 

. . n ,n 4Tsin 2 (7ra/A0 . , n , 9 1 . , 

11 Q 4QT sin 2 (ira/N) V ; 

The instability condition (|43[) and the spectrum of eigenvalues (|42"D give in- 
formation about dimensions of invariant manifolds around the uniformly os- 
cillating solution (a torus) characterized by \A°\. 
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Figure 1: The time sequences of three-site model, (a) trajectory constructed 
from the original map, (b) the naive RG map and (c) the regularized RG map. 
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Figure 2: (a) The phase space of the two-site model near the hyperbolic 
torus for \A°\ = \A C \ + 0.005. Here (0) corresponds to a uniform solution 
Aj = A exp(iQ\A°\ 2 n), and (1), (2) correspond to a™ and a^, respectively. 
The initial condition is as follows : ReA® = \A°\/V2 + 0.0002, ReA° = 
\A°\/V2 + 0.0001,1mA? = |A |/V2 + 0.0,imA° = \A°\/y/2 + 0.0. (b) The 
initial conditions are same as (a) when \A°\ = \A C \ + 0.5. 
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Figure 3: (a) The phase space of the three-site model near the hyperbolic 
torus for \A°\ = \A C \ +0.005. Here (0) corresponds to a uniform solution A r - = 



A exp(iQ\A°\ 2 n), and (1), • • •, (3) correspond to a™, 



respectively. The 



initial condition is as follows : ReA® = \A°\/y/2 + 0.0003, ReA° 2 = \A°\/^ + 
0.0002, ReA° 3 = \A°\/y/2 + 0.0001, ImA? = \A°\/V2 + 0.0, lmA° 2 = \A°\/y/2 + 
0.0,1mA!] = \A°\/V2 + 0.0. (b) The initial conditions are same as (a) when 
|A°| = \A C \ + 0.5. 




Figure 4: (a) The phase space of the four-site model near the hyperbolic 
torus for \A°\ = \A C \ + 0.005. Here (0) corresponds to a uniform solu- 
tion Aj = A exp(iQ\A°\ 2 n), and (1), • • •, (4) correspond to a", • • • , a\ , 
respectively. The initial condition is as follows : Re A? = \A°\/V2 + 
0.0004, ReAl = \A°\/V2 + 0.0003, ReA° 3 = \A°\/V2 + 0.0002, ReA° 4 = 
\A°\/V2 + 0.0001,1mA? = \A°\/V2 + 0.0,Iim4° = (A )/^ + 0.0, ImA° = 
1^°!/^ + 0.0,Iim4° = lA !/^ + 0.0,. (b) The initial conditions are same 
as (a) when \A°\ = \A C \ + 0.5. 
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Figure 5: The phase space of the three-site model near the hyperbolic 
torus for \A°\ = \A C \ + 0.1. Here (0) corresponds to a uniform solution A™ = 
A exp(iQ\A°\ 2 n), and (1), • • •, (3) correspond to a™, • • • , a^, respectively. The 
initial condition is as follows : ReA® = \A°\/\/2 + 0.003, ReA® = \ A°\/y/2 + 
0.002, ReAl = \A°\/V2 + 0.001, Iim4° = | A°\/V2 + 0.0, lmA° 2 = \A°\/y/2 + 
0.0,1mA!] = \A°\/V2 + 0.0,. 
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Figure 6: The time sequences of |a™| 2 of the three-site model near the 
hyperbolic torus for \A°\ = \A C \ + 0.1. Here (1), (2) and (0) correspond to 
| a" | 2 , \a 2 \ 2 and |ag| 2 , respectively. The orbit is characterized by the sequece 
of site-numbers 10210120 • • •. The initial condition is as follows : ReA® = 
\A°\/y/2 + 0.003, ReA° 2 = \A°\/V2 + 0.002, ReA° 3 = \A°\/y/2 + 0.001, ImA® = 
\A°\/V2 + 0.0,lmA° 2 = \A°\/V2 + 0.0,TmA% = \A°\/^/2 + 0.0. 
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Figure 7: The time sequences of |a™| 2 of the four-site model near the 
hyperbolic torus for \A°\ = \A C \ + 0.1. Here (1), (2), (3) and (0) cor- 
respond to | 1 2 and \a2\ 2 , respectively. The orbit is charac- 
terized by the sequece of site-numbers 112301 • ■ •. The initial condition is 
as follows : ReA\ = \A°\/y/2 + 0.004, ReA° 2 = \A°\/V2 + 0.003, ReA° 3 = 
\A°\/y/2 + 0.002, ReA° 3 = \A°\/y/2 + 0.001, ImA® = \A Q \/^2 + 0.0, lmA° 2 = 
\A°\/V2 + 0.0,lmA° 3 = |A°|/V2 + 0.0,ImA° = \A°\/y/2 + 0.0, 




Figure 8: (a) The Bernouilli sequences of the three-site model near the 
hyperbolic torus for \A°\ = \A C \ + 0.1. (b) same as (a) for four-site model. 
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